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Abstract 

The infrared behaviour of gluon and ghost propagators, ghost-ghion vertex and three-gluon ver- 
tex is investigated for both the covariant Landau and the non-covariant Coulomb gauge. Assuming 
infrared ghost dominance, we find a unique infrared exponent in the d = A Landau gauge, while 
in the d = 3 + 1 Coulomb gauge we find two different infrared exponents. We also show that a 
finite dressing of the ghost-gluon vertex has no influence on the infrared exponents. Finally, we 
determine the infrared behaviour of the three-gluon vertex analytically and calculate it numerically 
at the symmetric point in the Coulomb gauge. 
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I. INTRODUCTION 



In recent years there have been extensive non-perturbative studies of continuum Yang-Mills 
theory using Dyson-Schwinger equations in covariant Landau gauge, ref. 0, 0] and in the 
canonical quantisation approach in the Coulomb gauge, ref. [3] . In the latter Gaussian types 
of wave functionals have been used for a variational solution of the Yang-Mills Schrodinger 
equation for the vacuum 0, IM 0- Minimisation of the vacuum energy density gives rise 
to Dyson-Schwinger equations, which are very similar to the ones arising in the functional 
integral approach in covariant Landau gauge. Some of the relevant Green functions have 
also been calculated on the lattice in Landau gauge in refs. 0, 0, HI and in the Coulomb 
gauge in ref. The Green functions obtained by solving the Dyson-Schwinger equations 
to 1-loop order, are qualitatively very similar to the ones obtained in the lattice calculations, 
at least in the case of the Landau gauge. 

In the present paper we are interested in the infrared limit of the basic propagators and 
vertices arising in the canonical quantisation approach of Yang-Mills theory in the Coulomb 
gauge. With an appropriate choice of the vacuum wave functional, the corresponding gener- 
ating functional of the Green functions is structurally very similar to the one of the functional 
integral approach to Yang-Mills theory in Landau gauge, differing only in the number of rel- 
evant dimensions and in the precise form of the action. In both cases the infrared limit 
of the generating functional is governed by the ghost sector, which (up to the number of 
dimensions) is the same in both gauges. Therefore we can treat both approaches simulta- 
neously. Throughout the paper. Landau gauge will refer to the path integral quantisation 
approach whereas Coulomb gauge will refer to the canonical quantisation approach. 

Previously, the infrared limit of the gluon and ghost propagators in Landau and Coulomb 
gauges, were investigated in ref. 11, 3| and in ref. j^], respectively. In ref. |11|] two different 



solutions for the infrared exponents in Landau gauge were found, while in ref. using 
the angular approximation only one solution for the infrared exponents was found (in the 
canonical approach) in the Coulomb gauge. We also note, that in ref. the infrared 
behaviour of the ghost-gluon vertex in the Landau gauge has been studied. 

In this paper we perform a thorough infrared analysis of the ghost and gluon propagators 
as well as of the ghost-gluon and three-gluon vertices for both Landau and Coulomb gauge 
without resorting to the angular approximation. We discuss the validity of the previously 
known solutions for the infrared exponents of the propagators. The role of the ghost-gluon 
vertex in loop integrals is investigated and certain infrared limits of this vertex are explicitly 
calculated. Analytically, the infrared behaviour of the three-gluon vertex is determined at the 
symmetric point and for one vanishing external momentum. Numerically, we calculate the 
three-gluon vertex in the Coulomb gauge at the symmetric point over the whole momentum 
range. Given the Green functions in the infrared, we can confirm the self-consistency of these 
solutions. Also, the infrared fixed point of the running coupling for both the Coulomb and 
Landau gauge is calculated from the infrared behaviour of the Green functions considered. 

The paper is organised as follows: In the section IH] we present the infrared form of the 
generating functional of the Green function in Coulomb and Landau gauges. In section 
mil we perform the infrared analysis of the various Green functions: the ghost and gluon 



2 



propagators, the ghost-gluon vertex and finally the three-gluon vertex. Here, we also discuss 
the running coupling constant in the infrared. A short summary and our conclusions are 
given in section HVl Some mathematical details of our infrared analysis of the Green functions 
are presented in two appendices. 



II. THE GENERATING FUNCTIONAL FOR THE INFRARED 



The generating functional of the Green functions of Euclidean Yang- Mills theory defined by 
a Lagrange density Cym is given by 



Z[j\ / VAJ[A]eic^ 



d^xCYM{x)+ d^xfAx)Al{x) 



(1) 



where the functional integration is over gauge fields, which are restricted by a gauge condi- 
tion and J[Al\ denotes the corresponding Faddeev-Popov determinant. A common and for 
perturbation theory convenient gauge is the covariant Landau gauge 



9- A = , 

in which the Faddeev-Popov determinant is given by 

J[A] = dei{-d ■ D\A]) . 



(2) 



(3) 



Here, D[A\ = d + gA denotes the covariant derivative in the adjoint representation of the 
gauge group, i.e. A"''' = A'^f"''^^ with /"'"^ being the structure constants. 

The generating functional (0) in the Landau gauge is also the starting point for the 
derivation of the coupled set of Dyson-Schwinger equations, which to 1-loop level have been 
extensively investigated in recent years. For a review see ref. 

The Faddeev-Popov determinant represents the Jacobian of the transformation to the (curvi- 
linear) transverse "coordinates" A-*-, satisfying the gauge condition (j21). The appearance of 
J[y4] turns out to have crucial physical consequences. In order to pick a single gauge field 
out of the gauge orbit A^ = ^UDW , U G SU{Nc), it is necessary to choose configurations 
from within the (first) Gribov region Q or more precisely from the so-called fundamental 
modular region, which is a compact subset of the Gribov region and free of gauge copies. 

In the canonical quantisation approach to Yang- Mills theory, one uses the Weyl gauge Aq = 
to avoid the problems arising from a vanishing of the canonical momentum conjugate to Aq. 
Furthermore, Gauss' law, which here is a constraint on the wave functional to guarantee 
gauge invariance, is conveniently resolved in the Coulomb gauge defined by Eq. (j2)) for Aq = 
|3|. The Yang-Mills Schrodinger equation in the Coulomb gauge has been variationally solved 
in refs. ^] with the following ansatz for the vacuum wave functional 



^^[A]=AfJ-^[A]exp 



-- I d'[xy]A^{xMx,y)A^{y) 



(4) 
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where uj{x, y) is a variational kernel and A is a real parameter. The choice A > seems to be 
appropriate since it enhances field configurations near the Gribov horizon where J[A] = 0. 
It turned out that in 1-loop approximation to the Dyson-Schwinger equations arising from 
the minimisation of the energy density, the resulting infrared behaviour is independent of 
the value of A. Instead, it is the occurrence of the Faddeev-Popov determinant in the 
Hamiltonian that is crucial to the Dyson-Schwinger equation In order to investigate 
general features common to both the Coulomb and the Landau gauge, we now introduce a 
generating functional that facilitates the evaluation of expectation values of field operators 
in the Coulomb gauge^: 



zPij] = (vl/|exp 



- / d^x C{x) + / d^x]^{x)A'l{x) 



where 



C{x) 



d'^x'A1{x)uj{x,x')A1{x') 



(5) 



(6) 



Setting A = w.l.o.g., as mentioned above, a reinterpretation of the Faddeev-Popov determi- 
nant by means of Grassmann valued ghost fields becomes feasible and one can subsequently 
use common Dyson-Schwinger techniques to derive the equations for the Green functions. 
From Eqs. (0) and (0) it is evident that apart from approximations, the expectation values in 
Landau and Coulomb gauge follow from the same generating functionals by merely swapping 
the dimension (either d = 3 or = 4) and the respective actions. The infrared behaviour 
in the Coulomb gauge, in particular, is well described by a stochastic type of vacuum, as 
argued in ^^|. That is, setting \E'[74] = 1 will yield the correct infrared behaviour since it 
is dominated by the Faddeev-Popov determinant, i.e. by the curvature in orbit space. This 
circumstance, called "ghost dominance", corresponds to setting A = and £ = in Eq. ©. 
In the Landau gauge, ghost dominance was found as well i-e. setting Cym = will not 

affect the solution in the infrared. One is led to the conclusion that the infrared behaviour 
of the solutions of Dyson-Schwinger equations are the same in Coulomb and Landau gauge, 
if we consider d = 3 and d = 4, respectively. Therefore, calculating moments of the following 
generating functional that solely involves the Faddeev-Popov determinant,^ 



Z^'^^ [j, a, a]=X^ JvA j V[cc] exp 



d'^x{-cD[A\-dc + j-A + a-c + c-a) 



(7) 



one recovers the correct expectation values of spatial components of field operators in the 
Coulomb gauge by setting d = 3 and also the correct expectation values of field operators 
in the Landau gauge by choosing d = 4, as far as the infrared is concerned. 



^ The integration of the path integral should be restricted to the fundamental modular region (FMR) A C 
which is free of gauge copies. However, integrating over fl will yield the same expectation values [l^. 
Integration over il is understood in the following. 

^ Integrating in a compact region such as Q ensures convergence of the path integral |l3 • 
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FIG. 1: The (complete) Dyson-Schwinger equation for the ghost propagator, denoted by a dashed 
hne with a fuh blob. The curly line represents a connected gluon propagator and the vertex with 
an empty blob is a proper ghost-gluon vertex. 

III. INFRARED ANALYSIS OF GREEN FUNCTIONS 

The truncated set of Dyson-Schwinger equations (DSEs) that govern the Green functions of 
the theory defined by Eq. (|7j) can be solved analytically in the infrared limit. To do so, we 
make the ansatz that the propagators obey power laws in the infrared, and determine the 
values of the exponents. The investigation of the vertex functions, in particular the ghost- 
gluon and the three-gluon vertex then follows by investigating the corresponding Dyson- 
Schwinger equations. 

In the following it will be sensible to study Dyson-Schwinger equations in (i-dimensional 
Euclidean spacetime. One may then specify to the Green functions of the Yang-Mills vacuum 
in the Landau gauge by setting d = 4, or the Landau gauge high-temperature phase for c? = 3 
[Tfij l . In the Coulomb gauge, one can derive Dyson-Schwinger equations for equal-time Green 
functions which corresponds to the choice d = 3. 



A. Propagators 

In refs. a solution for the infrared behaviour of the propagators was previously 

obtained. We briefly review here the derivation of these results and critically analyse the 
approximations made. On the basis of our analysis we will have to discard one of the 
solutions found in ref. 
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The crucial dynamical properties of Yang-Mills theory are accessible via the computation 
of its two-point Green functions, i.e. the ghost propagator, 

D^^ip) := J d^'x {c\x)c\y)) e^^ ^--^) = ^^^'^^ , (8) 

as well as the gluon propagator, 

DtM ■■= I d'x {Al{x)Al{y)) e^^'^^"^) = = 5%.^)^ • (9) 

J P 

These expectation values as well as the set of DSEs that entangles them with higher n-point 
Green functions can be derived from the generating functional, see e.g. 0, 0]- Part of the 
information about the infrared behaviour of the ghost and the gluon propagator can be 
extracted from the ghost DSE. Using the bare ghost-gluon vertex^, it reads 

G-\p) = -^-gN, j d'^i (^l-ii-pf^Dzii)DG{i-p), (10) 



A discussion of a more general form of the vertex will follow further below. 
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where the integral is represented by the self-energy diagram shown in Fig. This integral 
bears an ultraviolet divergence for both three and four dimensions, due to the ultraviolet 
behaviour of the propagators. It can be conveniently subtracted by making use of the horizon 



condition 17 



limG'-^(p) = , (11) 

to find the finite expression 

G-\p) = G'\p) - G-\0) = gN, [ d'^i Z{i){l - {£ ■ pf)iDGii) ~ Ddi - p)) . (12) 



Aiming at the behaviour of G{p) for p — > 0, it is instructive to assume that below some 
intermediate momentum scale ^, the propagator dressing functions obey 

G{p) ^ G^-\p) = , Zip)^Z(-\p) = ^, p<C. (13) 

As done in refs. jH, 3|, the integral in Eq. ()12j) can be analytically evaluated if we naively 
replace the propagators in the integrand by the power laws given by Eq. (fT^ . even though 
the integration is over all space. This procedure is referred to as "infrared integral approxi- 
mation" from now on. It leads to the following representation of the integral equation ()12|) 
in the infrared 

(p^)-o = (p^)^-o-'^zAB^p^jG{aG, az) + $g(p') , (14) 

where Iq is a dimensionless number calculated in appendix A, see Eq. ()A14|) . and is the 

error to account for the infrared integral approximation. Neglecting this error, one obtains 
the "sum rule" 

d-A . . 

az + 2aG = , (15) 

along with 

AB^NJg{k) = 1 . (16) 

The infrared ghost exponent k, := ac can then be found by plugging Eq. (fT^ into the gluon 
DSE, as will be done below. 

Before, however, let us make an estimate for in order to understand which values of k 
are at all consistent with the sum rule ()15|) and when it is reasonable to make use of the 
infrared integral approximation. What happens when employing the replacements (|T3|l can 
be visualised by Fig. El For any generic self-energy type integral, 

d'^eD{e)D{i-p) , (17) 

the integration space can be divided into regions where the factors in the integrand yield the 
infrared power law and ones where they do not.^ Here, inside of a d- dimensional sphere of 



^ The following arguments can be applied as well to non-perturbative one-loop integrals of higher n-point 
functions. 
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FIG. 2: Sketch of the integration space of (|12|) for p < Inside the full circles, the propagators 
are given by their infrared behaviour. Adding and subtracting the infrared power laws outside the 
circles makes it possible to evaluate the integrals analytically. 



radius ^ around the origin, one has D{i) = D^"'\i), and inside a sphere displaced by p, one 
finds D{i—p) = D^'^^\i—p). The crucial point is that it is only inside of these spheres, where 
the integrand may become non-analytic, i.e. the ghost propagator has a pole at vanishing 
momentum, due to the horizon condition. As long as p <^ C,, we can always find a third 
sphere in the intersection of the two others, see Fig. |21 inside of which all possible poles of 
the integrand lie. If we then add and subtract D^^^\i) D^^'^\i — p) outside the third sphere, it 
becomes possible to represent the original integral / by a sum of integrals. Asymptotically, 
as p — s> and all three spheres intersect, this sum reads 

and we can define the first term plus the third term to be /q = / d'^i D^'^'^^i) D^^^\i ~ p) = 
/ — $, the infrared integral approximation of I. Jq is an integral which regards the integrand 
factors as the infrared power laws p3|) over all space and "captures" any singularities of the 
original integrand. The error of the approximation, $ = / — /q, integrates analytic functions 
and we can therefore expand it into a power series, 

oo 

$(p^) = ^a„(/r , p^O, (18) 

n=0 

which is finite at p = 0. The integral Jq, on the other hand, diverges as p — > for those values 
of the infrared exponents where Jq exists. This can be understood from the convergence 
criteria for two-point integrals given in appendix A. $ is then negligible for p ^ and one 
can set I = Iq. On the contrary, if p <^ ^ is not satisfied, the term $ does have a substantial 
contribution, as can be seen, e.g., in numerical calculations in section ITlI CI 
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FIG. 3: The (truncated) Dyson-Schwinger equation for the gluon propagator. In the infrared, 
other terms are neghgible. 

Returning to the ghost DSE ()12|). we note that the renormahsation plays an important role 
for the error <I>g of the infrared integral approximation. Within the subtraction of the UV 
divergence the term with n = cancels in a power series such as ()18|). We therefore find 
^cip^) = C^(p^) and infer that we can neglect this term for p — in the sum $g(p^) + (p^)"^ 
of Eq. (fT^ as long as 

< K = < 1 • (19) 

The lower bound is due to the horizon condition On condition of the above relation, 

the sum rule (fT^ is satisfied. For values of k with k > 1, the power series $g does have 
to be taken into account, and one arrives at a sum rule different from (jisp . However, those 
values are discarded since they do not allow for Fourier transformation of the propagators. 
Turning our attention to the gluon propagator DSE, we note that since the Faddeev-Popov 
determinant J'[A] dominates the infrared, only the ghost loop has to be included. This has 
been found in the Landau gauge |^ il^] as well as in the Coulomb gauge 01 for the equal- 
time gluon propagator. After contracting the gluon DSE with the transverse projector and 



taking the trace, we find [11 



Z~\p) = g'N,-^^^-^ I d'ef (l - (j3 ■ ir) Da{i)DGii-p) , (20) 

see Fig. El This integral is convergent in the ultraviolet. Employing the infrared integral 
approximation might introduce a spurious ultraviolet divergence, depending on the value of 
K. Of course, this divergence has to cancel with the error $2 of the approximation, since a 
choice of the infrared behaviour of the integrand will not affect the ultraviolet. Hence, there 
is no lower bound on k other than the horizon condition. The upper bound given by Eq. 
(fT^ guarantees convergence of the integral in the infrared. We then find 



ip'T' = ip')'^-''AB'NJziK) + ^z{p'), J9^0, (21) 

where /z(/t) is given in the appendix A, see Eq. ()A15|) . The error $2 is completely negligible 
since it approaches a finite constant in the infrared limit whereas the other terms diverge. 
This can be understood by noting az = (rf — 4)/2 — 2k < for k > 0. Thus, for p ^ 0, Eq. 
reproduces the sum rule (fT3|) and gives 

AB^NJziK) = 1 . (22) 

Along with Eq. (jl6j) . this leads to 

/g(«:) = Iz{^) , (23) 
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the conditional equation for k. For c? = 4, the Landau gauge case, only one unique solution 
lies in the range given by Eq. (0, k(^) ^ 0.595 [lil, Q- The solution k(^) = 1, claimed 
in jll|, could not be confirmed.^ In the case d = 3, applicable for the Coulomb gauge, we 
find two solutions, Kp"* ^ 0.398 and k^'^ = 1/2, in complete agreement with jll|. Only one 

of them, is found using the angular approximation |^. Numerical calculations in |^ 

(c) 

approximately approach the other value k\ . However, an improvement of the numerical 
methods shows that k'^'' is also a stable solution 3| • The latter would result in a Coulomb 



potential that rises strictly linearly. Which one of the solutions is energetically favoured is 
therefore an interesting issue that is yet to be investigated. 



B. Ghost-gluon vertex 

The infrared behaviour of the ghost-gluon vertex in the Landau gauge of SU{Nc) Yang-Mills 
theory for both d = 4 and d = 3 was investigated in Dyson- Schwinger studies 0, 3] and 



in lattice calculations for SU{2) and d = 4 j8|,|9|. It was found that its non-renormalisation 
which holds to all orders of perturbation theory jl^l, remains valid in the non-perturbative 
regime. An appropriate question to ask is, to what extent does a finite dressing of the ghost- 
gluon vertex influence the numerical solution for the infrared exponent k of the propagators? 
Leaving aside the colour structure which is assumed to be that of the bare vertex, i.e. /°*'^, we 
denote the proper reduced ghost-gluon vertex by r^(fc; where k is the outgoing gluon, q 
the outgoing ghost and p the incoming ghost momentum. Following ref . fl^] , a quite general 
ansatz for the ghost-gluon vertex is^ 

r„ft«.p) = .»Ec^.(^) (24) 

where the constraint li + nii + ni = 0, Vz, guarantees the independence of the renormalisation 
scale a, i.e. non-renormalisation of the vertex. It is readily shown that the sum rule 
fll5|) is not affected by a dressing of the ghost-gluon vertex such as ()24|) . since it turns into 

az + 2aG = + + mi + Ui) . (25) 

i 

Further investigations in 3| showed that the value for k, determined by Eq. (j23p . only 
slightly depends on the values of {L, m,;, rii}. 

Since neither the DSE studies [l^, I3| nor the lattice calculations jl, 0| show any infrared 
divergences, the dressing function of the ghost-gluon vertex must be some finite function. 
To investigate the consequences of a finite dressing function of the ghost-gluon vertex, let 
us assume, for simplicity, that it is given by a finite constant, 

V,{k-q,p) = CVf{q) , (26) 



^ The reason is that ^ = Um^^i T^-r4 

2 h.^^ Ig(k) 



7^ Vnnj^^4 ( hiiiK^i 777^ ) = 1- Furthermore, for k = 1 the term 

$G is not neghgible in Eq. (|14|) . 



,4 ■ \ Io{k)^ 

$G is not neghgible in Eq. (|14|) . 
^ GeneraUy, there is another component along the gluon momentum which, however, has no contribution 
when contracted with a gluon propagator. Therefore, it can be discarded. 
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where T^^\q) = igq^ is the bare ghost-gluon vertex. Then, the infrared analysis of the 
propagators can be performed in the same way as above. The ghost self-energy and the 
ghost loop are both multiplied by the constant C. In Eq. ((221); ^^is constant appears on 
both sides to one power and thus trivially cancels. Therefore, a constant dressing of the 
ghost-gluon vertex is completely irrelevant for the infrared behaviour of the propagators. 
The question that arises is if a non-constant dressing of the ghost-gluon vertex might result 
in a change for the determining equation ()2Hj) of n. The investigations in fisl] showed that 
after one iteration step of the ghost-gluon vertex DSE, the vertex remains approximately 
bare over the whole momentum range, i.e. C ~ 1. Also, the results in |13|] confirmed 
the well-known fact [i^ that for vanishing incoming ghost momentum p, the ghost- vertex 
becomes bare in the Landau gauge. ^ It is essential for the proof given in ref. [20'] that the 
gluon propagator is strictly transverse. It has been argued that the same holds true for the 



Coulomb gauge [21|, where the gluon propagator is transverse as well. If we discard the 
irrelevant component of the vertex along the gluon momentum /c, T ^ also becomes bare for 
vanishing outgoing ghost momentum q |l2,|l9|, i.e. 

\imT^{k-q,p) = \imT,{k-q,p) = Tf{q) . (27) 

However, the infrared limits of the ghost and gluon momenta are generally not interchange- 
able.^ In particular, zero gluon momentum yields a dressing that is different from one, 
although quite close to it, as we will see. The following relation shall redefine C, 

CrW(g)^lim(limr^(fc;g,p)) ^limflimr^(A;;g,p)') =Vf{q), q^O. (28) 

Does C 7^ 1 or a non-constant C affect Eq. (j^ ? To see this, it is not necessary to get 
involved in a numerical calculation but we can argue qualitatively instead. Consider any 
loop integral that involves the ghost-gluon vertex. Wherever it may appear in the loop 
diagram, the ghost-gluon vertex is always attached to ghost propagators. The integrand 
will be strongly enhanced for those loop momenta where the ghost propagator diverges, i.e. 
for p 0. Since the gluon propagator, on the other hand, is finite for all momenta, any 
infrared singularities in the integrand of the loop integral can actually be due to the ghost 



^ This agrees with the corresponding Slavnov- Taylor identity in the Landau gauge. 

® In this context, one might note that the infrared limit of any tensor integral is non-trivial. Given an 
integral 

we can construct a tensor basis from the external scales {p*-*'} and Lorentz invariant tensors. According to 
the Passarino-Veltman formalism, the above integral can then be expanded in this basis, which is nothing 
but solving a set of linear equations for the expansion coefficients in this basis. If one sets up a tensor 
expansion for finite {p^^^} and then tries to perform the infrared limit of a single external momentum, say 
pik) _^ ^jjg coefficient matrix becomes singular, and the tensor expansion is not well-defined. Instead, 
one can set p'-'"'^ — from the beginning (if the integral exists here) and construct the tensor basis 
spanning a vector space which is of a lower dimension than originally. The expansion coefficients are then 
well-defined. 
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FIG. 4: The (truncated) Dyson-Schwinger equation for the ghost-gluon vertex. 

propagator only. The ghost-gluon vertex does not introduce any additional singularities, 
since its dressing function is finite. For the value of the integral, the singularities in the 
integrand will give the dominant contribution. The only value of the dressing function of 
the ghost-gluon vertex that is relevant to the integral, is then the one where any of the ghost 
momenta vanish. According to Eq. the vertex is bare in these limits. We can therefore 
infer that in any loop integral the bare ghost-gluon vertex will yield the correct result. This 
circumstance can thus be traced back to the horizon condition and the transversality of the 
gluon propagator. 

Nevertheless, the constant C, defined by Eq. ()28|) is not entirely meaningless since the 
introduction of a running coupling, see section UlI Dl below, makes use of it. One can actually 
analytically calculate C by means of the DSE for the ghost-gluon vertex 3]; see Fig. |3J 



r,(fc;g,p) = rW(g) + rf^^)(A;;g,p) + rf^^)(fc;g,p) . (29) 
Here, F^f^*^^'' is a graph with two full ghost and one full gluon propagator in the loop, 

Ff^^)(fc;g,p) = dHY^^\-t)D^p{i- q)Yf,{q- t-v-t-k) 

DG{i + k)r^{k;e,i + k)DGii), (30) 

and Y^^^^^ has two gluon and one ghost propagator in the loop, but involves a proper 
reduced three-gluon vertex T^^p, 

N P 

Tf\k;q,p) = -Yj d''eT^^\q)DG{i-q)Tp{e + k;q-i,p) 

Dppii + k)r^,p{k; e, -i - k)D,^{i) . (31) 

Since F^(fc; q,p) exists in the limit — > jsl, Isl, (l3|, we set = in the integrands which 
greatly simplifies the tensor structure of Eqs. ()3()j) and (jH^- Furthermore, the proper ghost- 
gluon vertices that appear in the loop integrals are rendered bare, as discussed above. We 
then get 

rf ^^)(0; q, q) = zg'C^ J d''iLD.p{i - q)qp^,Dl{l) , (32) 



and 



rf ^^^0; g, q) = g'^ J d'^t q^qpD^,{e)Df,p{e)DG{i - g)r^.p(0; i, -£) . (33) 

Naively, we would expect from ghost dominance in the infrared that the contribution (jHHj) is 
subdominant since it incorporates only one and not two ghost propagators, like Using 
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a bare three-gluon vertex, we can calculate both integrals for g in the infrared integral 
approximation and indeed find that (j^Hj) becomes negligible. The dominant part of the two 
is then (jHH), and it gives (see appendix A) 

limrf«^)(0;g,g) 



which agrees exactly with the results of numerical calculations of (j^^ in this limit [l9|. 
Because of this agreement, we infer that the error introduced by the infrared integral ap- 
proximation, employed for the calculation of (j32j) and (j33|) . vanishes. 

If the dressed three-gluon vertex is included, see below, the graph (jH!^ has a substantial 
contribution to this limit of the ghost-gluon vertex. The calculation is then somewhat more 
involved, see appendix A, but one can extract the values for C in all cases at hand: 



1.108 


for d = 


4, 


K = 


«:(^) 


^ 0.595 


1.089 


for d = 


3, 


K, = 




^ 0.398 


1 


for d = 


3, 


K, = 


AC) 

1^2 


1 

~ 



(36) 



It is quite remarkable that in the Coulomb gauge with the solution Kg = 1/2, the two 
non-trivial graphs that appear in the DSE for the ghost-gluon vertex show an exact mutual 
cancellation in the infrared gluon limit, 

hm (rf ^^)(0; q, q) + if ^^)(0; g, g)) = . (37) 

Therefore, the interchangeability of limits is recovered in this case only and the ghost-gluon 
vertex becomes bare in all infrared limits. Note also that the results fl36|) are independent 
of Nc- The colour trace that occurs in the loop diagrams of Eq. yields a factor of Nc/2, 
see Eqs. (jHUjl and (jHT|) . but it cancels with the propagator coefficient term AB'^ = 1/{IgNc). 



N I 

Tf(q)'/C^ j dH e„g{q)t„„{l - q)Dl{e)Dz{t - q) (34) 

rrfa)c jf:^'...,.^-^:^r'^^:'\ , (35) 



f + 2/c) 



C. Three-gluon vertex 

As mentioned in section |nj it is only the Faddeev- Popov determinant that influences the 
infrared behaviour of Yang-Mills theory. The solution obtained for the infrared exponents of 
the propagators was found to be independent of the three-gluon vertex, in particular, since 
it does not contribute to the ghost self-energy term or the ghost loop contribution to the 
gluon self-energy. In the Coulomb gauge, we have seen that any of the vacua ^a[^] given 
by Eq. (0)) minimises the energy w.r.t. A, evaluated to one-loop order in the DSE (two-loop 
order in the diagrams for the energy). The question is how the three-gluon vertex changes 
in the infrared for different values of A without resorting to the one-loop approximation. 
The full three-gluon vertex is defined as 

^t>^y.z) = {m;i^)Aliy)A%z)\^) . (38) 
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FIG. 5: The (complete) DSE for the three-gluon vertex derived from the generating functional 
given in Eq. Q. 



In the particular case of A = 1/2 it is found that 



(39) 



by symmetry. Hence, the three-gluon vertex vanishes for ipi/2 |^. Now consider the case 
A 7^ 1/2. In [i^ it has been found that the Faddeev-Popov determinant can be written to 
one-loop order as 



J [A] = exp 



- / d''[xx']x^';{x,x')At{x)A){x') 



(40) 



with X being the curvature of the space of gauge orbits, i.e. the ghost loop contribution to 
the gluon self-energy. In this form, J' merely modifies the inverse gluon propagator u that 
occurs in Eq. (0) according to 



uj{x, x') — > n{x, x') := uj{x, x') - (2A - l)x(x, x') . (41) 

The subsequent determination of Vt by minimising the vacuum energy guarantees that Vt 
always is the same function, regardless of the A chosen. In particular, = for A = 1/2. 
Therefore, all expectation values are the same to one-loop order in the equation of motion, 
so the three-gluon vertex vanishes for any A. 

On the other hand, without the use of the one-loop approximation, A 7^ 1/2 will give a 
non-zero three-gluon vertex, in contrast to Eq. (jHUj) . as will be shown. Thus, the three- 
gluon vertex shows great sensitivity to the choice of the vacuum wave functional ip\i ^ 
behaviour not exhibited to one-loop order. Making the choice A = permits the standard 
representation of the Faddeev-Popov determinant by ghosts. In the following the infrared 
behaviour of the three-gluon vertex for this case will be explored following the treatment in 

23- 

The Dyson-Schwinger equation for the three-gluon vertex is derived in Appendix B and de- 
picted diagrammatically in Fig. Its complete form comprises a diagram with the unknown 
two-ghost-two-gluon vertex which is truncated here. The finite ghost-gluon vertex appears 
here in a loop integral and can therefore be set to its bare value throughout, according to the 
discussion given below Eq. in the last subsection. Assuming tree-level colour structure 
for all of the correlation functions and Fourier transforming the truncated DSE, one arrives 
at 

T,.,{puP2,P3) = N, [ d'eDG{k)DG{ps + k)DG{k-p2)T^;^\i)T^^\i-p2)T^^\i + Pi) (42) 
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pi «>(n) 

FIG. 6: Dressing function of the Coulomb gauge proper three-gluon vertex at the symmetric point. 
The dashed curve shows in contrast the perturbative case where the propagators in the loop are 
bare, i.e. k = 0. 



where the outgoing momenta obey the conservation law 

Pi + P2 + P3 = . (43) 

The vertex given by Eq. (ji^ is projected onto the tensor subspace spanned by the tensor 
components of the tree-level vertex. Due to Bose symmetry, the coefficient functions of these 
six components are all the same, but their signs alternate as the vertex without the colour 
structure is antisymmetric under gluon exchange. One finds 

^i,up{Pi,P2,P3) = - i{P2)t^SupF{pl,pl,pl) + i{p3)f,S^p F{pI,pI,pI) 

+ i{Pl)J^,pF{pl,pl,pl) -i{p3)uS^p F{pI,pI,pI) (44) 

- i{Pi)p^i^uF{pl,pl,pl) + i{p2)pSpu F{pI,pI,pI) . 

Equating Eq. (jl^ with (jl^ and contracting with these six tensors, yields a set of six linear 
equations for F, the solution of which reads 

_/Y r 

^^^''^''^'^ ^ io{pipi-{k-P2y) J ^'^^^w^g(p3+^)i^g(^-p2) 

{{pI+Pi ■ P2)(-2J2 -J4 + 3J5) + {pI+Pi ■ P2)(2Ji - 3J3 + Je)) (45) 

where 

J,: = {p,.i)f J,: = {p,.i)f J,:={p2-i){pi-i) 

J,:=plf J,: = {P2-i)' -HPi ■ P2)i' . ^ ^ 

The integral (jl^ depends only on the ghost propagator in this truncation, and despite 
the infrared enhancement of the latter it is convergent. The numerical calculation of the 
form factor F at the symmetric point, where pi = p\ = p\ ='■ p"^, shows a strong infrared 
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divergence, see Fig. IHl For the ghost propagator we used the numerical results of ref. |^ 
where k = 0.425 ^ k\ . A fit to the data in Fig. IHl yields an infrared power law such 
as F{p'^) ~ (p^)"^ ''^. If we employ the infrared integral approximation to calculate the 
symmetric point analytically, a power law behaviour can be extracted as well. By momentum 
scaling of the integration variable in Eq. (jl^j) . £ — > Xi, one finds Fip^) ~ (^p2^c(/2-2-3re_ 
Plugging in the value k = 0.425 obtained numerically in gives an infrared exponent 
of —1.775 which agrees (within errors) with the numerical result. This shows that the 
infrared integral approximation becomes exact as the external momenta vanish. Only in the 
ultraviolet there is a deviation from the infrared power law which can be clearly explained 
by the error of the approximation. At large momenta the vertex vanishes which complies 
with asymptotic freedom since in this approximation there is no tree-level vertex, due to 
Eq. ©. 

The above results for the infrared behaviour of the Coulomb gauge three-gluon vertex can be 
generalised to any value of k and any dimension d. Therefore, we can also make statements 
about the Landau gauge. Noting that d/2 — 2 — 3k = az — o^g, see Eq. (|T^ . we find 

Fip') ~ (47) 

With the analytical results for k in Coulomb [d = 3) as well as in Landau {d = 4) gauge, 
this yields 



F(/) 



^p2y.775 
1 



1 



for d 



for d 



K = K 



(C) 



(C) 



785 



for d = A , K, 



K 



0.398 

1 

2 

0.595 



(48) 



The Landau gauge result agrees exactly with ref. 

Another interesting kinematic point is the one where one of the gluon momenta, say pi, is 
set to zero while the others remain finite. Trying to calculate this point from Eq. ()44|) by 
setting pi = 0, the projections onto the tensor components fail as the determinant of the 
coefficient matrix that defines the tensor expansion vanishes in this case. It is advisory to 
impose pi = in the DSE (jl2D, 



r^,p(0,p, -p) = -tg'N.J d''ei^{i-p)JpDUi)DG{i-p) 



(49) 



One can then realize that this integral exists. It can be expanded into a tensor basis con- 
structed by the only scale p and Lorentz invariant tensors, i.e. {p^S^p, Pu^p^, PpSfj,u}- However, 
the only component that survives the transverse projections of the gluon propagators at- 
tached to the legs of F^^^p with finite momenta, is obviously Pp5yp. Thus, we can write 



^pypi^^P, -p) 



-ig^NcP^6, 



up 



id - l)p- 



-,Patp^{p) / d^i Upi^Dl{l)DG{l -p) + ... (50) 



where the ellipsis represents irrelevant tensor components which shall be discarded hence- 
forth. Using the infrared integral approximation, we find a finite expression, see in 



15 



appendix A, which makes the three-gluon vertex 

r,.p(0,P,-p) = -iB^N,p,5^, ^^^^[l_^^ , p^O. (51) 
The error can be ignored due to the infrared enhancement of (f^ . 

In view of the strong infrared divergence of the three-gluon vertex, one has to check for ghost 
dominance in the propagator DSEs, which simply states that one is to count the infrared 
exponents of the propagators in the loop However, the vertex function have to be taken 
into account, too [2]. The infrared power law of the three-gluon vertex ()47|) . expresses that 
the vertex dressing replaces the infrared exponent of a gluon by that of a ghost propagator, 
for any dimension d. The infrared hierarchy of terms in the gluon DSE remains untouched, 
since even with the dressing of three-gluon vertex, terms involving it remain subleading in 
the infrared. E.g., the gluon loop, which has an infrared exponent of (i/2 — 2 — 2az with a 
bare three-gluon vertex, attains an infrared power law with the exponent (i/2 — 2 — — ao 
if the vertex is dressed. Clearly, this term is still subleading w.r.t. the ghost loop which 
bears an infrared exponent of (i/2 — 2 — 20^^- 



D. Infrared fixed point of the running coupling 



A renormalisation group invariant that qualifies as a non-perturbative running coupling can 
be extracted from the ghost-gluon vertex and is given by 



"0 



Here, ao = S'o/^^ is the bare coupling constant and zl^'^{p' 
ghost field renormalisation functions, respectively. According to the definition of the renor- 
malisation function ^i(p^) of the ghost-gluon vertex in we find that Z^\0) = C, with 
C given by Eq. (jHUj) .^ In the Landau gauge {d = 4) this leads to an infrared fixed point 
ac = a{0), 



(52) 



z}J^{p'^) are the gluon and 



a 



8.915 



10.94 



(53) 



Up to the vertex correction C^, this value was found in l3|. Although the result holds 
for both Landau and the interpolating gauges j2l|, the Coulomb gauge limit reveals an 
infrared fixed point different from it, even if C = 1. By definition 2l|| one gets 

AAB'^C'^ 



a 



37T 



(54) 



With the values for the integral Ig{i^) at (i = 3 for the solution k^^^ 
correction = 1.187, this yields 



0.398, and the vertex 



a 



{C)\ 



11.99 



AC) 



14.21 



(55) 



^ Note that we have used the definition of C in Eq. H28() . With the infrared hmits not being interchangeable 
generally, an alternative renormalisation prescription of the ghost-gluon vertex would give C = 1. 
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For the other Coulomb gauge solution, where K2 =1/2 and the vertex correction vanishes, 
we find 




'2 



(C) — 



3iV, 



IGtt 



16.76 



(56) 



Among the interpolating gauges, the running coupling appears to have an infrared fixed 
point that changes discontinuously in the Coulomb gauge limit. 

A second possibility for a definition of a running coupling is given by the instantaneous 
Coulomb potential which is singular in the infrared jl5|. These two choices clearly disagree 
on the description of long-range interactions. It is known that the Coulomb string tension 
is an upper bound for the string tension extracted from the Wilson loop [2^ . 

IV. SUMMARY AND CONCLUSIONS 

We have studied the infrared limit of the ghost and gluon propagators as well as of the ghost- 
gluon and three-gluon vertices to 1-loop order in both Coulomb and Landau gauge assuming 
ghost (loop) dominance in the infrared. From the Dyson-Schwinger equations we have found 
that there is a unique infrared exponent in the Landau gauge, while there are two different 
exponents in the Coulomb gauge corresponding presumably to different minima of the energy 
density. It would be interesting to determine which one corresponds to the absolute minimum 
of the energy density. In the Coulomb gauge for one of the infrared exponents we found an 
exact cancellation between the two loop-diagrams of the DSE for the ghost-gluon vertex. 
This vertex is infrared finite in both Landau and Coulomb gauge. We have also shown, 
that a finite dressing of the ghost gluon vertex does not modify the infrared exponents. 
The three-gluon vertex was found to be infrared divergent with approximately the same 
infrared exponent in Coulomb and Landau gauges. Furthermore, the infrared divergence 
of the three-gluon vertex does not spoil the infrared dominance of the ghost loops over the 
gluon loops. We also calculated the three-gluon vertex numerically in Coulomb gauge over 
the whole momentum range for the symmetric point. The numerical result reproduces the 
infrared behaviour found analytically. For one vanishing gluon momentum, we determined 
the three-gluon vertex for any dimension. Furthermore, the three-gluon vertex turned out 
to be quite sensitive to the specific form of the Yang-Mills vacuum wave functional, and 
therefore a lattice calculation of this quantity would be of great interest. Finally, we have 
also calculated the infrared fixed point of the running coupling and found a larger value in 
Coulomb gauge than in Landau gauge. 
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APPENDIX A: TWO-POINT INTEGRALS 



Here we sketch the derivation of the integrals necessary to calculate the Feynman graphs in 
the infrared integral approximation. The two-point integrals, 

E^{a,P):=J ^^,^V((£_^)2)/3 ,«,/3eR,meN, (Al) 

which are encountered in the calculations can be shown to be homogeneous functions of the 
momentum q. By a scaling of the integration variable, i — > Xi, one readily finds that since 
the two-point integral can only depend on the scale, it should obey S^(q;,/9) ~ (?^)'* and 
the exponent of the power law can be determined to be 

d/2-a- P + m. (A2) 

After applying the usual trick of introducing Feynman parameters, 

1 

dxdy5{x + y - 1) , „ , , (A3) 



CfCf J ' {xC^+yC2Y+PB{a,(5)' 







where B{a, (3) is the Euler beta function, we can shift the integration variable, i ^ £ — yq. 
The integrand then depends on i"^ only and we can integrate it out. For m = 0, 1, 2, 3 we 
need the following standard integrals: 

d'e _ 1 r{n-d/2) fiy-'/' ^^^^ 



I 
I 



(£2 + A)" (47r)<^/2 r(n) \A 



(£2 + A)" 2 ''"(47r)'^/2 r(n) 



Integrals with an odd number of vectors £ in the numerator vanish by symmetry. For our 
purposes, wc have A = xy(f. The Feynman integrals can be straightforwardly computed 
using the identity 

1 



for the beta function. One then finds the results 



-o(«,/3)-^^ r(a)r(/3)r(d-a-/3) ^ ' ^^^^ 

^ R\ 1 r(d/2 - « + l)r(d/2 - /3)r(a + - d/2) . 

= r(a)r(/3)r(d-a-/3 + i) ^ ' ^^^^ 

^ rn, /^^ = 1 r(d/2 - g + 2)r(d/2 - /j)r(a + p- d/2) o 
^2ia,p; T{a)T{l3)T{d-a- P + 2) ' 

1 1 r(d/2 -a + l)r(rf/2 - /3 + l)r(a + /3 - ci/2 - 1) o^d/2-a-/3+2 ..q^ 

2(47r)'^/2 r(a)r(/3)r(c^-Q;-/? + 2) ^ ^ 
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„ 1 r(rf/2 -a + 3)r(c^/2 - [3)V[a + (3- d/2) . o.d/2-a-0+3 

^3l«,PJ T{a)T{(3)T{d-a- [3 + 3) ' 

3 1 T{d/2 -a + 2)T{d/2 - /? + l)T{a + ^-d/2-l) 2^./2-.-ff+3 ....^ 

2 (47r)'^/2 r(«)r(/5)r(rf - « - /3 + 2) ^ • ^ ^ 

The above formulae are valid for those values of a, /3, m only for which the integrals converge. 
At £ = g, a pole is integrable as long as (3 < d/2. On the other hand, the infrared convergence 
at £ = depends on m: 

J d/2 + m/2 for even m ^Aii^ 

^^\d/2 + m/2 + 1/2 for odd m ' ^ ^ 

One can relate this inequality to the requirement that the arguments of the first gamma 
functions in the numerators of ()A7|1 . ()A8|1 . ()A9|) and ()A10|) be positive. The ultraviolet 
convergence is contained in the third gamma function of the numerators. For convergence, 
the relation 

, . / d/2 + m/2 for even m , . 

" + ^ ^ \ c//2 + m/2-l/2 for odd m ^ ' 

has to be satisfied. Obviously, odd values of m, compared to even values, work "in favour" 
of convergence both in the infrared and in the ultraviolet, due to the angular integration. 
It is worthwhile examining two-point integrals that do not converge. Making the usual 
replacement 

t.q = l{^e + q'-{t-qf) (A13) 

in any of the convergent integrals Sm, with m > 0, one can express these in terms of a sum 
of So integrals, but generally encounters both IR and UV divergences. Curiously, if we use 
the regular result ()A7|1 for the divergent Sq integrals anyhow, the sum will yield the correct 
result given by the direct formulae ()A8|) . ()A9|) or ()A10|) . as can be checked. This indicates 
that divergent two-point integrals can be written as a regular part given by the formulae 
calculated here, plus the divergence which may cancel with another integral of the same 
kind. 

This circumstance is of great use for the ghost DSE, where the subtraction of G^^(O) removes 
the UV divergence and we can calculate /g('^) in Eq- ffT^ by 

IG{^^) = (P')"" (22(^ -2kA + k)/p^ - So(^ - 2k - 1, 1 + z^; 

V / reg. 

A-{d-i) r(f-K)r(-K)r(i + K) ^^^^^ 



(47r)'^/2+V2 r(f - 2K)r(l + f + K 



The same result was found in [ll| where subtractions of divergences were circumvented. 
The integral that occurs in the gluon DSE ipiH) is essentially Iz{,i^) defined by Eq. (j2H), and 
can be calculated to give 

Iz{k) = (p2)-«_l_(H^(^^i + ^)_H2(l + «:,l + «:)/p2) 

1 r(f-Kfr(i-f + 2K) ^^^^^ 



2(47r)'^/2 r((i - 2 k) r(i + k) 



2 
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Note that although both integrals in ()A15|) have an infrared divergence at i = q, the sum is 
regular for those values of k in ()19p . as can be seen by shifting i ^ i — q. In the ultraviolet, 
flA15|) only converges as long as k > d/4 — 1/2. However, such a divergence will necessarily 
cancel with the error as discussed in section 1111 Al Regardless, the solutions for both 
the Coulomb and the Landau gauge yield k > d/i — 1/2. 

For the calculation of the infrared limit of the ghost-gluon vertex, we encounter the integral 
in dnni)- With 



L,{qMi -q) = - -^0^, (A16) 



we find, using Eq. 



limrf«^)(0;g,g) = Vf{q)Cg'^AB- + 2/., ^ - 2k) - -3(2 + 2«:, ^ - 2^)/^^) 

= rf(g)cif (A17) 

where 

_1 d — \ ( K^Q^ 

1 - (47r)<i/2 d (1 + 2K)r(rf/2) ■ ^ ^ 

Plugging in the value ()A15|1 for Iz{n), leads directly to Eq. (jH^j) . 

Before one can calculate F/t (0; q, q), the three-gluon vertex is needed. According to Eqs. 
(ISUj) and (f^. we can find the integral Is to yield 

^ - ^(Si(l + 2K,l + K)-S3(2 + 2K,l + «:)/p2) 



1 F(f-2K)F(f-K)F(2-f + 3fi:) 



2(47r)'^/2(c/ - 1) F(rf-3K)F(l + /t)F(2 + 2K) ' 
Using this expression for the three-gluon vertex in Eq. ()33p. we find that 

limFf^^)(0;g,g) = -^g^A'B'hq,^ J d'i i ■ q {l - {p ■ If) Dlmdi - p){p 



(A19) 



2\-{ac-az) 



-igq^ i ^ (Hi(rf/2 - k, 1 + k) - ^3(^/2 - k + 1, 1 + K)/q^) 
z 

-Tf{^)Y-j^ (A20) 



where 

1 ci- 1 



(47r)'^/2 (rf - 2K)V{d/2 + 1) 



(A21) 
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Altogether, the infrared gluon hmit of the ghost-gluon vertex, defined by Eq. ()28|) yields 



^"^ + 2^/. 2 11 



= 2 {4'' id-l)d (r(l + r(^)) T{d - 3«:)r(^ - «:)r(-K)^r(i + «:)r(l + 2«:) 

r(2 + 2/t) - r(^ - 2/€)=^r(i + ^ + /€)2r(2 - ^ + 3«:) 

{(d-l) {d~2K) T{d -3 k) T{-k) 1(1 + 2 fi;)^ 

(4^^'^r(i + ^)r(^ - K)r(-«:)r(^ + «:) + v^r(^ - 2«:)r(i + ^ + '^)) } (A22) 

As can be checked, this leads to the numerical values of C given by Eq. ()36|) for the various 
solutions of K. 



APPENDIX B: THE DSE FOR THE THREE-GLUON VERTEX 

The Dyson-Schwinger equation is derived from the generating functional Z of the theory. 



Z[j,a,a]= j V[Acc]e^-p - j d'^xC + j rf'^x ^A)^ + a'^c'^ + cV) 
where C is given by 

C{x)= [ d'^x'A''^{x)uj{x,x')A''^{x') - [ d'^x'c''{x){-d- D[A]y\x,x')c\x') 



(Bl) 



(B2) 



see Eq. ©, with A = 0. 

To derive the DSE we observe that 







V[ccA] 



- / AujA+ j c{-d-D[A])c+ J (jA+ac+ca) 



(B3) 



as the integral can be turned into a surface integral over the Gribov horizon where the 
Faddeev-Popov determinant vanishes. We perform the derivative and replace emerging 
fields by derivative operators with respect to their sources in order to recover the generating 
functional Z. It is replaced according to Z = e^, thus introducing the generating functional 
W[j, a, a] of the connected Green functions: 



-2 / d X uj(u, x) 



5j 



Here we use 



^(x) J ^ ^6a^{x)^ ' 

6i-d-D[A]){y,z 



6A-{x) 



and 



(5^ (5^ 

— = derivative acts from the left; — = derivative acts from the right . 
OCT oa 



(B4) 



(B5) 



(B6) 
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The derivative in Eq. ()B4|) corresponds to a first gluon of tlie tliree-gluon vertex. For tlie 
otlier two gluons we perform two further derivatives, 



and 



(B7) 



on Eq. ()B4|) . Setting j = cr = a = in Eq. ()B4|) and in the equations obtained after each of 
the derivatives in Eq. ()B7|1 yields the DSE for the connected one-, two- and three-point gluon 
Green functions respectively. Plugging the one- and two-point DSE into the three-point DSE 
we get 



2 / d'^x u{u,x)- 



+ I d'^[xx']{Tl;')''\u;x,x') 



5^W 



Sj^{w)6j^{v)6a^{x')6a''{x) 



m 



where J denotes all the sources collectively. After the application of 



d'^z u~\y,z)uj{z,x) = S'^iy -x) (B9) 
we can identify the gluon propag ator D^l on the RHS, cf. Eq. (37) in Introducing 



WZ{x,y,z) :- 



WT'iw.x.y.z) :-- 



Sj;{x)6jt{y)6j-{z) 
6^W 



S3?,{w)6f^{x)6cx^{y)6a^{z) 



(BIO) 



j=o 



provides as an intermediate result for the DSE for the connected three-gluon Green function: 
Wt^(x,2/,^) = [ d''[uvw]Df^{z,u)iTyrfiu;v,w)W^y^{x,y,w,v) (Bll) 



To derive a DSE for the proper three-gluon vertex we decompose the connected Green 
functions into the proper ones by using the generating functional F of the proper Green 
functions which is the functional Legendre transform of W defined by 

T[A, c, c] := -W[j, a, a] + / d'^x (A;(x)j;(x) + a^xyix) + c»(x)a"(x)) (B12) 



where the sources on the right side are chosen to fulfil 



Al{x) 



Sa^ix) 



c^i^x) 



5W 

da'^ix) 



(B13) 



Therefrom we derive the relations 



5T 



5A-{x) 



50" ix) 



5c'' ix) 



(B14) 
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and the inversion relation 
d z- 



5^T 



'Sjt{y)Sj-{z)6A;{z)6A-{x) J 6jt{y) SA;{z] 
which is the starting point for the decomposition. We apply 

S 



(B15) 



(B16) 



to Eq. ()B15jl , use the chain rule and further apply Eq. ()B13| IB14| IB15jl . With the definitions 



jjab 



and 



J=0 



6A-{x)6Ab{y)6A';,{z) 
the decomposition of the connected three-gluon Green function reads 

W^t;i^,y,z) = - [ d''[uvw]V^i{u,v,w)D^^{u,x)Dt{v,y)Di'^^{w,z) , 



j=o 



(B17) 



(B18) 



which simply means cutting off the external propagators. This is different in the decompo- 
sition of the two-ghost-two-gluon vertex. The starting point is an inversion relation similar 
to Eq. ()B15|) that also is derived from the Legendre transformation: 



d'^z- 



6^W 



6^T 



,Ja\y) 5c%z) 



d^z 



S''^5'^{x - y) 



(B19) 



5a''{x)5(j''{z) 5c''{z)5&{y) J ~ dc^iz) 5a''{x) 

To actually perform the decomposition it is useful to view the second derivatives as matrices 
where the colour, Lorentz and coordinate indices are matrix indices and the integration is 
the matrix multiplication. So we write Eq. ()B19|) as 



5a5a 



5c5c 



-1 



From jj.{AA ^) = 0, where A is a matrix, we obtain 



5A-^ 



A-' . 



6t " 5t 

We use this in taking the derivative of Eq. ()B20|) with respect to the gluon source: 



5^W 



6j^{x)5(75a 



6c6c 



6^r 



6jl^{x)6c6c 



6c6c 



(B20) 



(B21) 



(B22) 



We apply Eq. ()B20|) and Eq. ()B13|) to this and take a further derivative with respect to 
the gluon source. Together with Eq. ()B22l IB18|) and making frequent use of the techniques 
just developed we obtain the decomposition of the connected two-gluon-two-ghost Green 
function. With the definitions 



Dh\x,y):- 



6^W 



6a''{x)6a^{y) 



J=0 



6Al{x)6c\y)6c''{z) 



6^r 



^ ' 6A-{w)6A>;,{x)6c^{y)6c^{z) 



(B23) 



j=o 
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and switching back to the index notation it reads 

d'[pqTstu]D<^^{s, z)D%{x, t)D;l{w, p)D^^{y, q)Dj!{r, u)T':^{p, g, r)r;^^(t, u, s) 
+ I d''[pqrstu]D^a\^, z)D^J{y, t)Df,{w, p)Dt\{x, q)Df^{u, r)T%{p, q, r)rf^{u, t, s) 
+ [ d'ipqrstujD^'is, z)D;1{w, t)D'J:^{x, p)D^^{y, q)D^J{r, u)T':^{p, g, r)Tl^^{t, u, s) 
d'istuvm'is, z)D%{x, t)D;i{w, u)D^a\y. ^)r^«r(«, t, v, s) . 



(B24) 

Plugging Eq. ()B24|) and ()B18|) into the DSE for the connected functions, Eq. ()B11|) . and 
solving for the proper three-gluon vertex yields the DSE for it. It contains, however, a 
term that seems to be one-particle reducible. This comes about as we have not yet taken 
into account the DSE for the lower proper correlation functions. The two-gluon DSE is the 
equation obtained after the first of the two derivatives in Eq. ()B7|) . Treating it the same 
way as the three-gluon DSE till the point of one-particle irreducibility gives 

d''[viV2WiW2]{T''^''r'{x;vuV2)D^Ji^2,wi)D^^\w2,vi)T^^^^ = , (B25) 

which actually is found to be a part of the improper diagram in the three-gluon DSE that 
consequently vanishes. 

We have thus arrived at the Dyson-Schwinger equation for the three-gluon vertex: 

r;%{x,y,z) = 

■T';f\x,v^,V2)T'Jfiy,w^,W2) (B26) 
+ / d''[uiU2V,V2] {Tl''Y'{z- u,,U2)D{f{v2, ui)D'J{u2, Vi)^;^^^ {x, y, Vi, V2) 
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